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Abstract 

The (001) surfaces of cubic SiC were investigated with ab-initio molecular 
dynamics simulations. We show that C-terminated surfaces can have different 
c(2x2) and p(2xl) reconstructions, depending on preparation conditions and 
thermal treatment, and we suggest experimental probes to identify the various 
reconstructed geometries. Furthermore we show that Si-terminated surfaces 
exhibit a p(2xl) reconstruction at T = 0, whereas above room temperature 
they oscillate between a dimer row and an ideal geometry below 500 K, and 
sample several patterns including a c(4x2) above 500 K. 
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Silicon Carbide is an attractive material for high temperature micro- and optoelectronic 
devices |]1| because of its wide band gap, high thermal conductivity, high hardness and 
chemical inertness. Having a small lattice mismatch with GaN, SiC has also emerged as a 
promising substrate for the growth of nitride-based devices 

In the last decade, a notable effort has been devoted to the characterization of SiC sur- 
faces, since SiC films for use in technological applications are prepared by epitaxial growth. 
In particular, the (001) surfaces of the cubic polytype /9-SiC have been studied with a variety 
of experimental [3-11] and theoretical [12-15] techniques. At the end of the eighties, it was 
established that these surfaces are terminated by only one species and a clear assignment of 
different LEED patterns to either C- or Si-terminated surfaces was given Nevertheless 
the reconstruction of the Si-terminated surface is still the subject of debate, in view of the 



partial disagreement between recent [10| and older [Big experiments. Furthermore, a long- 



standing controversy [H,|,12-14[ about the structure of C-terminated surfaces was solved 



only recently [|Tl|]. At present, important problems such as the effect of preparation condi- 
tions on the surface structure and the thermal stability of the (001) surfaces are yet unsolved. 
These are key issues for the understanding of any growth process on SiC substrates, and 
represent as well fundamental problems in the physics of compound semiconductor surfaces. 

In order to address these issues, we performed a series of ab-initio molecular dynamics 
(MD) simulations []l6l of several /?-SiC (001) surfaces at finite temperature. We studied 
various reconstruction paths of the C- and Si-terminated surfaces, modeling different prepa- 
ration conditions, and we identified experimental probes which could discriminate between 
the various reconstruction patterns. Our computations constitute the first ab-initio analysis 
of the effect of preparation conditions and thermal treatmtent on the structure of (001) SiC 
surfaces. 

Our calculations were carried out within the local density functional approximation 
(LDA). We used slabs periodically repeated in the (x,y) plane (see Fig. 1), terminated 
on each side by the same atomic species and followed by a vacuum region of ^ 8.5 A. In 
most MD simulations we used 11 layers, each containing 8 atoms. Structural optimizations 



at r = were carried out both with 11 and 19 layers (152 atoms) IjT^. The lateral di- 
mensions of the supercell were 8.60 x 8.60 A, corresponding to the theoretical equilibrium 
lattice constant of bulk SiC (4.30 A). The experimental value is 4.36 A. We note that the 
use of a symmetric slab is essential to avoid spurious charge transfers and electric fields, 
which could be present in H-terminated non symmetric slabs or in slabs terminated 



by different atomic species on opposite sides ||1J]. The interaction between ionic cores and 



valence electrons was described by fully non-local pseudopotentials |T8| with s and p non 



locality for Si and s-only non locality for C. Single particle orbitals and charge densities were 
expanded in plane waves with kinetic energy cutoffs (-Ecut) of 40 and 160 Ry, respectively. 
We considered Bloch functions at the T point of the supercell surface Brillouin zone (SBZ); 
this corresponds to including the points F, M, S, and F, J' and {J+J'/2) of the c(2x2) and 
p(2xl) primitive cells (001) SBZ, respectively. 

C-terminated surfaces - Experimentally, C-terminated surfaces are obtained either by 
ethylene decomposition on Si-terminated substrates or by Si sublimation be- 

tween 1200 and 1500 K. The two techniques lead to samples characterized by different 
diffraction data ^ and photoemission spectra We devised two series of computations, 
to investigate the two preparation methods. 

The first series of calculations concerns Si sublimation experiments. Since Si removal 
is expected to yield an ideally C-terminated surface (C-SiC(OOl)), which then reconstructs, 
we studied the spontaneous reconstruction of an ideally terminated C-SiC(OOl) and the 
thermal stability of its reconstructed phases. Starting from the ideal geometry vibrating 
at room temperature, we performed microcanonical MD simulations. Depending on the 
initial conditions on ionic positions and velocities, the ideal geometry transforms either into 
a c(2x2) staggered dimer (SD) pattern, or into a p(2xl) dimer row (DR) geometry (see 
Fig. 1). At T = 0, the total energy of the DR is lower than that of the SD geometry by 



0.44 eV/surface dimer |T9|. We note that the C dimers of the first layer have a similar 



bond length in the two reconstructions (1.36 and 1.37 A in the DR and SD, respectively). 



very close to that of the dimers on the diamond (001) surface pO|. In both geometries. 



the Si-C distance between the first and second layers is sensibly longer than its bulk value 
(~5%). However, in the SD reconstruction only the first C layer shows important deviations 
from the ideal geometry, whereas in the DR reconstruction three layers take part in the 
reconstruction, with both the second and third layer having a large buckhng ( 0.2 and 0.37 
A, respectively). 

To investigate whether a C-SiC(OOl) surface can exhibit an SD and/or DR reconstruction 
in the temperature regime where Si desorption occurs, we studied the thermal stability of 
the SD and DR geometries. In our calculations, the SD geometry is stable against thermal 
fluctuations up to 1200 K. At this temperature the distance between the flrst and second 
layer becomes considerably larger than at T = 0. Eventually some of the Si-C bonds 
break, yielding C dimers bonded only at one end to two Si atoms of the second layer. 
After some oscillations, these dimers rotate by about 90 degrees, and are flnally bonded 
to Si atoms at both ends with one bond per C atom. Si atoms in the second layer form 
dimers. We call such a configuration a bridge (B) pattern (see Fig. 1). We also obtained 
a B reconstruction by deposition of C dimers on a Si-SiC(OOl) surface, as we will discuss 
below. Our simulations clearly show that an SD geometry can transform to a B geometry 
at high temperature, whereas no tendency to form DR configurations was observed within 
the spanned temperature range. 

When heating the DR reconstruction up to 1800 K, we did not observe any topological 
change of the surface in spite of a considerable elongation and occasional breaking of the 
Si-C bonds between the uppermost layers. The C dimers mostly oscillate in the surface 
plane. Dimerization of some Si atoms in the second layer was observed above 1200 K. In 
the DR geometry, 90 degree dimer rotations that could give rise to B configurations are less 
favored than in the SD geometry; only a collective motion of all the dimers could allow such 
a structural transition. 

In summary, depending on initial conditions, an ideally terminated C-SiC(OOl) surface 
at finite temperature can spontaneously reconstruct into two different geometries, which 
are both stable up to the temperature range of Si sublimation experiments. Therefore the 
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reconstruction resulting from Si desorption is expected to be determined by the kinetics 
of the subhmation process and the precise temperature at which Si removal occurs, rather 
than by the energetics of the surface at T = 0. Furthermore we find that the coexistence 
of various reconstruction patterns, including DR, SD and B geometries, is possible in the 
temperature range of sublimation experiments |21| ]. 



The photoemission data of Parril et al. [^] for Si sublimated surfaces did not reveal 
the presence of any Si-Si bonds, which would instead be present in a B geometry. This is 
consistent with our simulations at T < 1200 K. Our results are also consistent with the data 
of Semond et al. [Q, who obtained rather inhomogeneous surfaces by thermal desorption. 
However there is so far no experimental evidence of p(2xl) reconstructions resulting from 
Si sublimation. This could be explained either by the kinetics of the Si desorption process, 
or by entropic effects at high temperature, which may favor c(2x2) reconstructions. 

In order to make contact with experiments preparing C-SiC(OOl) surfaces by ethylene 
deposition on Si-SiC(OOl) substrates, we performed a second series of computations. We 
arranged C dimers in a staggered pattern on an ideally Si-terminated surface and optimized 
the geometry of the system. The dimers were given two inequivalent (perpendicular) ori- 
entations on the top and on the bottom surfaces. After annealing, a B reconstruction was 
obtained on both faces: the top C dimers rearranged so as to induce a marked dimerization 
of the second layer Si atoms. The bottom C dimers turned by 90 degrees to yield even- 
tually the same geometrical pattern as on the top surface. We obtain a remarkably small 
(1.23 A) C-C distance in the C dimers, indicating triple bonds, and a bond length of the 
Si dimers (2.37 A) similar to that in bulk Si (2.35 A). Contrary to what is found in the 
SD and DR reconstructions, the Si-C distance between the uppermost layers decreases by 
about 3% with respect to its bulk value. As in the DR geometry, three layers take part in 
the reconstruction. 

The calculated total energy of the B reconstruction is lower than that of the SD by 
0.18 eV/dimer and higher than that of the DR by 0.2 eV/dimer. These results agree well 
with those of Sabisch et al. [l^, although our energy differences are larger than in their 



calculations, and are consistent with the findings of Ref. . 

Our simulations, showing that C2 molecules deposited on a Si-SiC(OOl) surface can yield 
a stable B reconstruction, confirm the experimental results of Long et al. |Tl|], who recently 



demonstrated that C2H4 deposition on Si-SiC(OOl) lead to a B geometry. It would be very 
interesting to explore whether other hydrocarbons (containing e.g. two-fold coordinated 
C atoms) might lead to the so far unobserved DR reconstruction. Powers et al. |^ also 
reported that ethylene deposition on Si-SiC(OOl) leads to a B geometry. However they 
suggested that Si desorption should yield as well a B reconstruction, in spite of difficulties 
encountered in fitting either an SD or a B geometry to their data for Si sublimated surfaces. 
These difficulties can be understood in view of our results, by a possible coexistence of SD 
and B configurations in surfaces prepared by Si desorption. 

We now turn to the discussion of the electronic properties of the C-SiC(OOl) surfaces 
and in particular of surface band structures (SBS). In the B reconstruction, the valence 
band top (VBT), located at F, has a clear bulk character, whereas in other parts of the 
SBZ the highest occupied state has a surface character. This surface state is a combination 
of (J bonding orbitals on the Si dimers and of tt^ bonding orbitals on the C dimers. The 
conduction band bottom (CBB) (also at F) has again bulk character. Above the CBB, 
we found two surface-bulk resonances. The lowest in energy has large components on the 
C dimers, given by antibonding combinations of orbitals, which are orthogonal to the 
dimer orientation (tt* state). At the SBZ boundaries, a surface state with antibonding 
components on C dimers was observed (p* state). These results are in very good agreement 



with the spectra of Long et al. which show a prominent Pz peak with surface character 



and a less intense vr^. state at higher energies [22|. As pointed out in Ref. the existence 



of a Tlx state is a clear signature of the B reconstruction. Furthermore Long et al. found 
evidence of states coming from Si-Si bonds from Si 2p CLS measurements. In our calculation 
the B structure is semiconducting with a gap of ~ 1.4 eV. 

The SBS of the SD and DR geometries are simpler than that of the B reconstruction, 
due to the absence of Si-Si bonds. In both surfaces the highest occupied (HO) and lowest 
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unoccupied (LU) states at the F point have bulk character. Close to the VBT and CBB, 



three surface states were observed: a tt^, a tt* and a tt* state, localized on the C dimers p3 



Our SBS results agree very well with those reported by Sabisch et al. |[15[| . 

In order to predict the results of STM experiments for the C-SiC(OOl) surfaces, we report 
in Fig. 2 the square moduli of the HO and LU surface states for the three geometries, on 
a plane parallel to the surface at a distance of 1.5 A. There is an unequivocal difference 
between the images of the p(2xl) and the c(2x2) surfaces. Drawing a clear distinction 
between the B and SD images is more delicate. When rotated by 90 degrees, the images 
of the vr* and vr^ states of the B surface may not be clearly distinguishable from those of 
the TT* and tt* of the SD surface. The key difference between the two sets of images lies in 
the state occupation: the n states perpendicular to the surface are occupied (empty) and 
those parallel to the surface are empty (occupied) in the B (SD) geometry. Fig. 2 suggests 
that STM experiments could be effective probes to identify the various reconstructions, and 
possibly reveal the coexistence of domains with different reconstructions on the same surface, 
as was recently observed in diamond (001) p4[. 

We further analyzed the differences between the reconstructions of C-SiC(OOl), by com- 
puting Si-2p and C-ls core level shifts with respect to bulk SiC. We adopted the procedure of 
Ref. and ionized atoms of the two uppermost layers ||2^. The values obtained for Si-2p 
CLS in the SD and DR reconstructions are similar (-0.12 eV and -0.09 eV), but different in 
magnitude and sign from that calculated for the B reconstruction (0.32 eV). Indeed, Si-Si 
bonds only exist in the B geometry. C-ls CLS of the atoms of the first layer are instead 
very similar in the three structures (-0.44 and -0.48 and -0.46 eV for the SD, DR and B 
geometries, respectively), although only in the B reconstruction C dimers are triply bonded. 



Our calculated C-ls CLS for the B geometry is smaller than the experimental value |11 



-1.1 eV), the error being close to that found in a variety of carbon systems ||2^. Whereas 



the LDA is believed to give a good quantitative description of CLS involving Si atoms [25 



only a qualitative agreement between theory and experiment is obtained for CLS of several 
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C systems For the B geometry we also computed the CLS of C atoms of inner layers. 
These CLS have magnitudes similar to those of the uppermost C atoms but have opposite 
signs. We therefore expect the presence of broad peaks in C-ls core absorption data. 

Si-terminated surfaces - In our calculations for Si-SiC(OOl), the ground state geometry at 
T = is a p(2xl) dimer row reconstruction, which lies a few meV/atom lower than the ideal 
surface geometry. The Si dimer bond length (dsi) is 2.58 A, larger than the value (2.31 A) 
fitted to LEED data 0. We investigated the thermal stability of the reconstructed surface by 
performing constant temperature MD simulations in the range 200 K < T < 800 K. At T ~ 
300 K the surface dimer length varies between ~ 2.3 and 2.9 A, i.e. close to the value of the 
ideal structure (3.03 A). At higher temperatures (300 K < T < 500 K), the surface oscillates 
between a p(2 x 1) dimer row reconstruction and an ideal surface geometry. The large thermal 
fluctuations of dsi could be responsible for the difficulties encountered in achieving a good fit 
of LEED data for Si-SiC(OOl) surfaces PJ7|JT0[| . At T > 500K, we observed a tendency to form 
a c(4x2) reconstruction, where neighboring dimers have opposite buckling [Q. The ability 
of the Si-SiC(OOl) surface to sample several reconstructions in a relatively small temperature 
range could explain the differences between recent |]TD[ and older LEED experiments [0,0. 
Finally, we evaluated Si-2p CLS and found a value (0.85 eV) which is in very good agreement 
with experiment (0.8-1.0 eV) |T0[| . 

In summary we investigated the effect of preparation conditions on C-terminated 
SiC (001) surfaces, as well as the thermal stability of both C- and Si-terminated surfaces. 
Computations of band structures and C-ls and Si-2p core level shifts of the stable C-SiC(OOl) 
geometries point at STM images and Si-2p CLS as effective probes to identify the different 
reconstructions. 

We thank A. Pasquarello and M. Sabisch for useful discussions. The calculations were 
performed on the NEC SX3/4 of CSCS in Manno, Switzerland. Work supported by the 
Swiss NSF Grant No. 20-39528.93 (GG) and the itahan CNR (AC) through a FNS-CNR 
cooperation. 
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FIGURES 

FIG. 1. Top view of the dimer row (DR), bridge (B) and staggered dimer (SD) reconstructed 
geometries of the C-SiC(OOl) surface. Three layers are shown. Black and white spheres indicate C 
and Si atoms. Axes are oriented as in Ref. [11]. In the three geometries dimers are not buckled. 

FIG. 2. Square moduli of the highest occupied (lower panel) and lowest unoccupied (upper 
panel) surface states in a plane parallel to the surface, at a distance of 1.5 A, for the DR (left, tt* 
and TT* states), the B (center, vr^ and tt* states) and the SD (right, tt* and vr* states) reconstructions 
of the C-SiC(OOl) surface (see Fig. 1). White bars indicate C dimers. In the B reconstruction, at 
1.5 A from the surface the electronic states are dominated by C dimers components and Si dimers 
components are not visible. 
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